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Abstract. We consider a family of surfaces of revolution ranging between a disc and 
a hemisphere, that is spherical caps. For this family, we study the spectral density 
in the ray limit and arrive at a trace formula with geodesic polygons describing the 
spectral fluctuations. When the caps approach the hemisphere the spectrum becomes 
equally spaced and highly degenerate whereas the derived trace formula breaks down. 
We discuss its divergence and also derive a different trace formula for this hemispherical 
case. We next turn to perturbative corrections in the wave number where the work 
in the literature is done for either flat domains or curved without boundaries. In the 
present case, we calculate the leading correction explicitly and incorporate it into the 
semiclassical expression for the fluctuating part of the spectral density. To the best of 
our knowledge, this is the first calculation of such perturbative corrections in the case 
of curvature and boundary. 



1. Introduction 

Waves can be approximated by rays by means of asymptotic expansions. For example 
in the design of optical lenses one always resorts to the equations describing the rays 
of light rather than the Maxwell equations for the electromagnetic waves. Similarly for 
the interpretation of earthquakes, one replaces elastic waves with curved rays. 

The parameter of such an expansion for a classical wave phenomenon is the wave 
length divided by a typical geometrical scale of the systems. In optics, this leads from 
electromagnetic waves to rays of light. In quantum mechanics, however, the asymptotic 
parameter is h compared to a typical action of the system, that is, it involves a natural 
constant. Nevertheless on the formal mathematical level, the asymptotic expansions 
have the same character and many semiclassical results have their analogue for classical 
wave equations in terms of similar ray interpretations. Likewise, there are trends in 
the opposite direction, from classical waves to quantum mechanics, such as imaging by 
electron waves in e.g. photodetachment microscopy pQ. For this article we stress that 
the semiclassical discussion of a free quantum mechanical particle on a curved surface is 
similar to that of the ray theory for a classical membrane. Thus, in curved spaces, the 
effect of geometry on quantum mechanics is that wave packets follow geodesic motion 
in the semiclassical limit in the absence of potentials j2]. Geodesies are curves which 
are stationary with respect to the length under variations. 
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In this contribution we consider a one-parameter family of spherical caps. We show 
that detailed information about the spectra is obtained from periodic geodesic orbits. 
In particular, we explain the drastic change in the spectrum involving a clustering effect 
when approaching the half-sphere by varying the opening angle of the cap. Besides 
the method of stationary phase, we shall use the asymptotic technique attributed to 
Jeffreys, Wentzel, Kramers and Brillouin and denote it JWKB. We shall discuss how to 
take into account the effect of the boundary with respect to the first JWKB-correction 
using the scattering approach of [3]. 

Our interest in spherical caps is motivated by the ongoing experimental efforts for 
elastic shell caps by Ellegaard and his group. 




R 



e 



Figure 1. The cap geometry: opening angle 9q and curvature radius R 

The article is organized as follows: to aid a classical interpretation we recall various 
facts about geodesic polygons in section |2j Then the spectral problem is introduced in 
section [3] with a discussion of the asymptotics of the associated Legendre polynomials. 
Using a scattering approach to quantization, this allows a derivation of a trace formula 
for the density of states, see section |4j Then follows a discussion of the pure hemisphere 
in section [5] which can be solved exactly along with a comparison in section [6] with the 
previous general case. We include JWKB-corrections in the trace formula in section [7j 
We end with a short discussion in section [8j Finally, calculational details are given in 
appendices. 

2. Classical quantities: geodesic polygons 



We fix the notation in section 2A_ and discuss the conservation of angular momentum 
during the geodesic flow in section 2.2 For reference later on, we compile a list 
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of geometric quantities for geodesic polygons in section 2J3 give a condition for the 
inclination of a segment in a geodesic polygon in section 2.4 and discuss the concept 
of anholonomy present in systems with continuous symmetries such as U(l) with a 



calculation for the case of the cap in section 2.5 



2.1. Notation 

We shall consider a family of surfaces of revolution. Individual members of this family 
are spherical caps and are parameterized by a radius of curvature R and an opening 
angle #0 calculated from the axis of symmetry taken to be the z-axis, see figure [TJ We 
consider caps ranging from almost a disc to a hemisphere and we restrict ourselves to 
caps with opening angles less than or equal ir/2. The singular cases corresponding to 
#0 = 0, n are not considered. 



2.2. Conservation of angular momentum 

If a surface possesses symmetries these influence the dynamics of the geodesic flow. For 
example, for surfaces of revolution Clairaut [5] found from the geodesic equations: 

rcos0 = Cst. , (1) 

where r is the radial distance from the rotation axis and is the angle of the tangent 
vector of the geodesic with a given latitude, see [5]. When Q is multiplied with the 
momentum p, the corresponding physical interpretation is that the angular momentum 
along the rotation axis is conserved. Here, one uses the fact that the momentum is 
conserved in the absence of potentials [2]. 

In the presence of a boundary, the geodesic flow is broken but it shall be continued 
by allowing geodesies to follow the law of reflection. Therefore, for the cap geometry 
the periodic orbits will be geodesic polygons composed of great circles. In the present 
case, the reflections are operations which preserve angular momentum. 

Likewise, from a spectral point of view, the angular momentum is a good quantum 
number. 



2.3. Geometrical orbit quantities 

We label each geodesic polygon with n edges and winding number M [f] with a pair of 
numbers (n, M). The pentagram orbit for example has (n,M) = (5,2), see figure [2] 
Consider a single edge in such a polygon. Let its vertices be given in direction by the 
two unit vectors ei and e 2 with 

ei = sin 6> x + cos 9q z (2) 

and 



e 2 = sin 6*o cos 5<p x + sin 8 sin 5<p y + cos 8q z , 
\ i.e. how many times the orbit winds around the north pole 



(3) 
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Figure 2. The pentagram orbit for opening angle 8q = ir/3. 



where the azimuthal angle traversed is denoted by 



n 

So the angle from the origin between the two points defining the segment 



(4) 



X = Z(ei,e 2 ) (5) 
can be found by e.g. 

cos x = e i ■ e 2 = cos v + sin tf cos • (o) 

n 

Hence the length of this polygonal edge is 

Al = R X (7) 
corresponding to one n'th of the total length of the orbit. 

2.4- Orbit inclinations from conservation of angular momentum 

With respect to the origin a particle moving on the segment with momentum p has an 
angular momentum L. Classically |L| and L z are conserved. The latter is also conserved 
under reflections from the boundary. The direction becomes 

L = T 2 -i (8) 

|ei x e 2 | 

and the magnitude is 

\L\=Rp. (9) 

Using ^ and ([3]), we find after some calculation the projection of the angular momentum 
on the z-axis: 



L z = pR cosip (10) 
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with 



sin # cos ^ , 
cos i> = n . (11) 

W#o + sin 2 # cos 2 ^ 



In the following, we shall interpret (11) as a condition fulfilled by an orbit's inclination 
(measured relative the polar axis) when given winding number M, number of bounces 
n and the opening angle of the cap 9 . 

2.5. Anholonomy 

We consider the geodesic flow as a dynamical system. We shall record the position of 
the trajectory at the boundary by a corresponding azimuthal angle and think of it 
as a phase. For an initial incidence angle this phase will change by constant discrete 
increments Sep and only if this phase is given by (J4J) the orbit is closed. We now enquire 
how this phase changes as the direction of the orbit changes. This direction we choose to 
be controlled by L z , the classical variable conjugate to the azimuthal variable <fi. Thus, 
if a general azimuthal increment 5(f) changes, L z changes via @, @ and After 
some calculation this yields 

d5(j> = 2 (cos 2 g + sin 2 g cos 2 (f)) 3 / 2 
dL z ~ pR sin 6 cos 2 6 sin (f) 
the phase lag [TElllE], also called the anholonomy. 

3. Semiclassics using asymptotics of Legendre polynomials 

We have finished the classical considerations and will next study spherical caps using 
wave mechanics. First, we shall discuss the wave problem of a spherical surface 
and the corresponding symmetry reduced problem given by Legendre's equation in 



section 3.1 We find the ray-limit of this reduced wave equation using the J WKB- method 



in section 3.2 This will lead to asymptotics for the associated Legendre polynomials. 



3.1. Wave problem of scalar cap 

The simplest model of cap vibrations corresponds to that of a quantized particle confined 
in the cap region: 

(A + A; 2 )^ = (13) 
with A the curved Laplacian, k the wave number and specified boundary conditions such 



as Dirichlet ^ = or Neumann d n ^> = at 6 = 6$. Equation (13) is also considered to 
describe curved drums. 



When solving (13) by separation of variables using 



m = u(ey m4> (14) 
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Legendre's equation 

„ 9n d 2 u du , m 2 \ , , 

( 1 - I )^- 2l S + ( i(( + 1 '-T^)" = ° < 15 > 

arises with the variable x = cos 9 and the parameter 

/(/ + 1) = (kR) 2 = k 2 (16) 

related to the spectrum. The regular solutions at 9 = are given by the associated 
Legendre polynomials -P/™(x) whereas the irregular ones by Q™(x). Consequently the 
Dirichlet solutions are the solutions of 

PTM = o (17) 

with respect to I and where xq = cos^o, see figure [3] showing log(|P z m (o:)|) in the case 



m = 20 and 9q = tt/3. When I < m the dips show an integer spacing. There as m is 
an integer PJ n {x) = for all x and integer I < m, since in general at integer l,m the 
Pf 1 is proportional to an m'th derivative of an Tth order polynomial. These "roots" 
are indeed present in figure [3] up to I = m — 1 = 19 but do not correspond to genuine 
eigenmodes as the wave function is globally zero. Around I ~ m there is a gap after 
which physical states appear spaced apart typically with a non-integer value for high I. 
The reason for this fixed spacing is due to the following asymptotics when I 3> m [6] : 

PTiccs 9) = r( ^ + 1} (vr sin 9)^ 2 cos [{l + 00 - \ + m|) + 0(1^) . (18) 



To further understand the roots of condition (17) in the interface region where I and m 



are comparable we turn to a more detailed asymptotics of the Legendre polynomials. 
3.2. JWKB of Legendre polynomials 

Standard tables of functions often discuss the asymptotics of Pf 1 ^) i n the extreme cases 
where either I > m oi m > I. However, as mentioned we need the more interesting 
intermediate physical situation with both I and m large and their ratio well defined 
following [7] (more explanatory details can be found in [HI El )• 

Thus, consider large angular momenta / and from that define a small parameter 

e = (/(/ + l))- 1 / 2 ^^ (19) 

and demand that the ^-component of the angular momentum m to be of the same 
asymptotic order as I. Hence, define the projected angular momentum on the z-axis 
respective x-y-plane 

m = L z and l ± = ^1(1 + 1) -m 2 . (20) 

The classical picture is the precession of the angular momentum around the z-axis at 



the angle ip, where we defined the corresponding cosine in (11). Thus for semiclassics 
we shall associate a classical angle ip to the quantum mechanical angular momentum by 

cosip = em and sin if) = e l± = a . (21) 




To simplify the following analysis, switch from the original wave function u to 

w{x) = Vl -x 2 u(x) . (22) 
Then Legendre's equation MM becomes 



9 d 2 w a 2 — x 2 + e 2 . , 

£ V + ir3^» = ° < 23 > 

free from the first order derivative term. This wave equation contains besides the kinetic 
term w"(x) also potentials: in the notation of [H] the leading potential is Q with: 

n 2 — t 2 

«•) = -IT^f < 24) 

and a sub dominant potential Q2 proportional to e 2 

1 

;i ~x i / 

To leading order in the JWKB-method this subdominant potential is dropped although 
we will also consider the effect of Q2 later on. 
First there is a momentum 

1-x 2 



Q^) = -7^-^- 2 - (25) 



p(x) = y^W) = ^ . (26) 



Second there is a turning point when this momentum vanishes, i.e. at x — ±a. We only 
consider < x < x < a. 

4. Semiclassical scattering quantization 



In the previous section we discussed the ray-limit of the symmetry reduced problem of 
a spherical cap in terms the asymptotics of the associated Legendre polynomials. In 
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this section, this information is used to derive the corresponding full spectrum with a 
ray- limit involving geodesic polygons. First, we state a simple asymptotic condition for 



an eigenmode in a form similar to that of Bohr-Sommerfeld quantization in section 4.1 



Then we reformulate the resonance condition as one arising from scattering in section 4.2 



and check this condition numerically in section 4.3 We next progress from this condition 



for individual eigenfrequencies belonging to a given irreducible representation to the full 



spectral density in section 4.4 We find, that the derived asymptotic density of states 



diverges in the limit of a hemisphere in section 4.5 



4-1. Resonance condition in the ray limit 

Following steps similar to the usual for Bohr-Sommerfeld quantization [TT] discussed in 



detail in Appendix A the Dirichlet condition \1/ = gives the following JWKB-condition 



for an eigenmode: 

\ 

e 

for n G Z and an action integral Iq given as 



\ + \ = nn (27) 



pit) dt = [ — dt = acos ( — 1 — Vl — a 2 atan ( ^ . =^— \ ■ (28) 

This condition follows from finding the roots of the asymptotic form of the wave function 

u ~ Y lm (0, <P = 0) « ^(a 2 - x 2 )- 1 ^ cos - ^ (29) 

taken proportional to a spherical harmonic. We derive this leading asymptotic form in 
Appendix A| as well as perturbative corrections for use later on in section [7j 



4-2. Scattering formulation 

We aim at deriving the spectral density for the spherical cap with the spectrum given 



by the condition (27) using ideas from scattering theory [3|. 

What is the scattering approach to quantization? This method exploits the 
connection between an exterior scattering problem and an associated interior resonator 
problem. Thus, if we force an obstacle from the exterior at an interior eigenfrequency 
the scattered wave experiences no phase shift; we shall see this explicitly in our case. 
The method is often referred to as inside-outside duality. We will apply this method to 
find the spectral density for the interior cap vibration problem using information from 
corresponding scattering data. We proceed to our particular geometry. There, formal 
scattering for curved billiards on e.g. the sphere has already been discussed in [T2]. For 



concreteness, we discuss the explicit form of the scattering states in detail in |Appendix 



Thus in terms of scattering theory, we shall say that a condition for an eigenmode 
is that the phase Q m : 

Q m =h 0+ 7 L (30) 
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has to equal zero modulo 2ir 

9 m = fmod27r) 



(31) 



First, (31) is formally equivalent to (27); second, the discussion in Appendix B shows 



that there is an underlying set of scattering states with precisely the scattering phase 
m . So we define a scattering matrix diagonal in the azimuthal quantum number m 

§ = (S m ) me z and S m = exp(-z 6 m ) . (32) 
Then an eigenmode occurs when the corresponding scattering problem is transparent 



[3J. The choice of sign in (32) agrees with treatments of the flat case, i.e. the disc [3] [T3]. 



4-3. Numerical test 

As an example, there is an eigenmode for I 
approximation for 9 = 7r/3, i.e. P; m (cos6'o) 



99.6428945787050 and m = 70 to good 
; 0. For these values of / and m we can 



find e and a from (19) and (21). The corresponding semiclassical phase evaluates to 



m / (27r) = 4.0014, close to an integer. 
4-4- Spectral density 

The distribution function for the eigenfrequencies is captured by the density of states. 
For convenience, instead of frequencies the spectral parameter in this section is the 
dimensionless wave number k = kR. For several spectral problems the density can be 
approximated by a decomposition consisting of a smooth p and oscillating part p. The 
smooth part is the most studied [H] whereas less is known about the oscillatory part. 
For our system we turn the attention to the latter. 

In the scattering formulation [3] the fluctuating spectral density becomes 



P{k) 



j oo 

— Im > 

7T ^ 



n=l 



1 d_ 
n dn 



Tr (S n ) 



— Im > 

7T ^ 



n=l 



1 d_ 

n dn 



Tr (S* r 



(33) 



The terms for the n'th power in (33) correspond to orbits bouncing n times as we shall 
see in the following. 



4-4-1- Poisson summation and n'th trace As the scattering matrix (32) is diagonal the 
trace becomes 



Tr §" = J2( S mV = $>xp(-m9 m ) ■ 

m m 

In the general case p4| ) is done by Poisson summation 

/oo 
dm exp(? (-n6 m - 2nmM)) 
■oo 



(34) 



(35) 



and subsequently approximated by stationary phase. 
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4-4 -2. The stationary point corresponds to geodesic polygons The saddle point 
condition for the trace (35) becomes 



— 71- 



dm 



2nM . 



We calculate the derivative of (|30|) using (|28|) 

-e atan 



dlo 
dm 



■s/o 



Xn 



a 2 xq 



-e atan 



a/ sin 2 ip — cos 2 9q 
cos ip cos 6*o 



and find the condition 

7fM 



atan 



a/ sin 2 ip — cos 2 #0 



// \ cos ?/> cos $0 

or by elementary trigonometric manipulations 

sin On cos — 

u n 



and 



COS?/' 



sin?/' 



cos 2 # + sin 2 9 cos 2 ^ 

COS #0 

cos 2 # + sin 2 #o cos 2 ^ 



(36) 



(37) 



(3* 



(39) 



(40) 



These conditions are precisely those of classical geodesic polygons winding M times 



around the north pole and hitting the boundary n times, see (11). For an example of a 



classical polygon orbit see figure [2] and figure [5] showing the pentagram orbit. 

4-4-S. Action at stationary point We next discuss the constant term in the stationary 



phase approximation of (35). 



The square of the denominator of (40) is rewritten using (|6 

2ttM\ 



In -In 2 nM 2 n ■ 2 n 1 + COS V „ 

cos Uq + sin t) cos = cos Vq + sin u — - 



n 



1 + cos x 2 X 

= cos — 

2 2 



with x being the angle of a single geodesic segment ([5]). Then consider 



acos 



cos 6*o 
sin?/' 



irM , x 

COS?/' = cosip 

n 2 n 



(41) 



(42) 



by (28), (40) and M. We proceed to calculate -nQ m - 2nmM in (35) by (30): the 



2-nmM cancels with the remaining term from (42) by (21). Thus 

7T 



2-nmM 



X 

-n n— 

e 2 



(43) 



For the trace formula (33) eventually the conjugate scattering phase is used 



corresponding to minus (43). By (|7j), nx/e = nxkR corresponds to the wave number 
times the total length of the geodesic polygon, i.e. the classical action of the geodesic 
polygon. At this point we identify the classical angular momentum Rp with multiples 
of the dimensionless wave mechanical, i.e. \Jl{l + 1), with pR replaced by kR. 
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4-4-4- Second variation From (37) the second derivative becomes 
d 2 I 2 x 



(cos #0 + sin #0 cos 



2 ttM\3/2 
n ' 



dm 2 



■ I n 



- n cos 2 8 sin 6q sin 

when evaluated at the saddle point. 

We express this Hessian using the classical anholonomy: 

d 2 I e dd(j) 



nM 
n 



(44) 



dm 2 



2 dL 7 



(45) 



when comparing (44) with (12). Such identifications are customary in the context of 



trace formulae for systems with continuous symmetries [T5l H~6] . 



As for the interpretation of the point of stationary phase in section 4.4.3 



we 



identify the classical angular momentum in (45) with the dimensionless wave mechanical 



momentum kR. Alternatively, we could have chosen a dimensionless measure of orbit 



anholonomy with a particle of unit classical momentum: d 5 (f)/d (cos ip), see (12). 



4-4-5. Trace formula Collecting the previous results leads to the following oscillatory 
density of states in the variable k = kR: 



p(k) 



EE 



dL 2 



X (—1)™ cos ( n\K — n 



71 



(46) 



M=l n=l 

We summarize the various factors in this formula. In particular x being the angle of a 
segment from the origin, n\n = k ■ nR\ is the classical phase proportional to the length 
of the geodesic polygon. The sign (— l) n is the phase shift from the Dirichlet boundary 
condition. Furthermore the derivative dL z /dA(p is interpreted as the accumulated phase 
lag for a closed orbit with n segments A0 = n5cj) using (12). 



The result (46) for the density of states of the spherical cap is consistent with the 



general £7(l)-symmetry reduced trace formulae discussed in [I5|[T6]. Similar agreements 
are found in the flat case for a disc, both quantum and elastic [31 fTTj . 



Equation (46) represents the result to leading order in k. Later in section 7 we 



include the first perturbative correction. 

The classical interpretation of the spectral density also holds more generally: much 
earlier in systems without symmetries Gutzwiller connected the fluctuations to classical 
periodic orbits [IB] using the method of path integrals. 

For a opening angles 8$ larger than tt/2 one can always associate the spectrum to 
that of a cap with opening angle tt — 6q as 







prr 



-X) 



x 



(47) 



with x = cos 6*o, already remarked in [12]. In particular, the geodesic polygons of a small 
cap also governs the spectral fluctuations of its complementary cap. 

Numerical checks: By including a few orbits and many repeats it is possible to 
build up a distribution with sharp peaks, see figure |4j This technique is discussed by 
e.g. [19] for the flat case of a disc. We have performed this simple test as a check for 
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Figure 4. Spherical cap: Oscillatory density of states for opening angle 6q = ir/3 



both Dirichlet and Neumann. In all cases at 6 = tt/3 the peaks fell approximately at 
the positions of the exact eigenfrequencies when including sufficiently many short orbits 
and repeats. 

4-5. Analysis of orbits as the opening angle changes 



We now consider the trace formula (46) for general opening angles and discuss what 



happens as the opening angle changes. As the opening angle is varied the set of periodic 
orbits changes accordingly. Figure [5] shows how the pentagram orbit gradually converges 
from a pentagram in the plane towards an orbit going twice around the equator. Likewise 
the semiclassical weights associated to each orbit change, see figure [5j The calculations 
were done using the first 14 orbits having winding numbers (2, 1), (3, 1), ... , (9, 2), (9, 4) 
and up to 80 repeats. 

Thus, the prefactor for each orbit family in the trace formula is (dL z / dAcf)) 1 / 2 . 

7r, A<f) = 2tt being the angular increments for the diameter orbit: 

dL z n ^ 



Using (12) where 8q 

{kR)- 1 ' 2 ■ 



( dLz \ 



de- 1 ' 2 + o{dd z ' 2 ) . 



(48) 



On the other hand non-diameter orbits have weights going rapidly to zero as a function 
of the deficit in the opening angle: 

dL 2 N 1/2 



(kR) 



-1/2 



(OLA 



O(d0) 



(49) 



V ^ ^ $ J non— diam 

In conclusion, in the limit of the hemispherical cap the diameter orbit controls the 
spectral density. This will also show up in section |5j 
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Figure 5. The pentagram orbit and the density of states for fixed radius (R 
and varying opening angles £ {5.7°, 26.5°, 47.2°, 68.0°, 88.7°}. 
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5. The hemisphere 

In the hemispherical case, the spectrum can be found exactly: we consider the spectrum 



in terms of individual levels in section [5TT| respective in the form of a spectral density in 
section I 



5.1. Exact spectrum of hemisphere 

For Dirichlet or Neumann conditions the eigenfunction in the polar variable either is 
odd or even. The parity of the Legendre polynomial is governed by 

Thus Dirichlet or Neumann conditions hold when either 

I — m = odd or I — m = even . (51) 

with degeneracies vd = I respective vn = I + 1. Thus the spectra will display bunching 
of eigenmodes at integer I respective approximate half integer k. 

5.2. Hemisphere density of states 




Figure 6. The hemisphere: orbits with the length of the perimeter. 

As we know the spectrum, the spectral density can likewise be written exactly. 
Denoting the degeneracy as v the density becomes 

oo 
n=0 
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En=o l6 ( l - n ) (Dirichlet) 

EZo( l + 1 ) S ( l - n ) (Neumann). 
Poisson summation then gives the density in the Dirichlet case 

oo . 

PD (i)= Yl i^ m + L -m 

N=-oo 

oo 



I + 21 Y cos(27riV/) 



N=l 



respective in the Neumann case 

Pn(1)= £ (l + l)^ Nl + l -±±5(l) 

N=-oo 

00 

= (/ + 1) + 2(1 + 1) V cos(27riV/) + - 5(0 . 

Af=l 

Switching to the spectral parameter k = fci? = / + 1/2 + 0(/ _1 ) for which koIk = 
(I + 1/2) all holds exactly gives 

« « " ^ + 2 - cos ^vriV - 0) (52) 

and 

Pjv(k)«« + ^ + 2 £ cos ( 27riV ( K +0) • ( 53 ) 

The first two terms above gives the smooth contribution to the density of states. 
Likewise the last terms correspond to the oscillatory part. 

The smooth counting function N is the integrated density of states and becomes 

^7 k2R2 kR 

N D/N = — T — (54) 

the first term being the available phase space volume divided by (27r) 2 whereas the latter 
is the boundary correction well known from flat Helmholtz resonators =|=fc/(47r) L with L 
being the length [TU [20]. Their sum corresponds to the spectral density of a full sphere. 
In particular, the sum of the boundary corrections add up to zero corresponding to no 
boundary. 

The fluctuations are governed by the phase 2nNn = k ■ 2ttRN of the iV'th repeat 
of an orbit with length of the circumference 2nR, see figure [6j 



6. Comparing the two trace formulae 



In the previous section 4.4.5 and section |5.2 we presented two trace formulae for the 
oscillatory density of states for a spherical cap. In this section we compare these two 
formulae. 
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Orbits: The first trace formula (46) derived for a general cap opening angle contains 



a countable infinite number of orbits whereas the second (52) only has a single orbit 
with the length of the perimeter. 

Weight versus wave number. In the hemisphere case this orbit has a weight 
proportional the spectral parameter and hence is of the order of the smooth part of 
the density. However, in the non-hemisphere case families of orbits have weights only 
proportional to the square root of the spectral parameter. This is the conventional result 
for families of orbits in two-dimensional systems. 



Weight versus opening angle: Inspection of the first trace formula (46) in section 4.5 



revealed that non-diameter orbits are assigned weights that vanish when the opening 
angle converges to tt/2, i.e. that of a hemispherical cap. On the other hand, the 
diameter orbit has a weight which goes to infinity in this limit. In conclusion the first 
trace formula is a non-uniform asymptotic result valid for < 8q < ir/2. On the other 
hand the second trace formula is only valid at 6 = tt/2. 



1000 



500 



-500 




Figure 7. Semiclassical density of states. Red/full line: 8q = 89.08° and blue/dashed 
line: 9 = 89.95°. 



Near hemisphere: When the opening angle approaches ir/2 (46) has been checked 
to reproduce the exact peak positions, k « 1 + 1/2 with i 6 N: a little bit away e.g. 
#o = 89.08° is not sufficient whereas at e.g. 9 = 89.95° the predicted half integers 
occur, see figure |6j This comparison was made at relatively high wave numbers with 
k = kR m 45. 

7. Leading correction 

The previous calculations were all done to leading order in the expansion parameter 
proportional to the inverse wave number e = (kR)^ 1 . We shall now ask what is the 
correction to these results in this expansion parameter? There are two points of interest: 



the individual sub-spectra for each m (section 7.1) and the total spectrum on the level 



of the trace formula (section 7.2) 



Besides the general question of how to calculate JWKB-corrections, we comment 
on why that could be of interest: one benefit of JWKB-corrections is that these may 
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increase the validity of trace formulae further down in the spectrum. Ultimately JWKB- 
corrections are believed to be particularly important in dimensions larger than two [3]. 
General results for corrections are known in the case without boundaries j2TJ E2] whereas 
only little is known with boundaries (planar Dirichlet /Neumann case [23]). By studying 
the spherical cap we will get direct access to a simple situation for a curved manifold 
with boundary. In the present contribution we shall obtain a non-trivial result already 
for Dirichlet conditions. 



7.1. Correction to leading asymptotics of associated Legendre polynomial 

We are interested in the correction to the action / when going from the boundary 



x to the turning point a for the equation (23). The discussion of JWKB-corrections 



when the orbit is between two turning points [TT] is our starting point. We outline 



the corresponding theory below with details of the calculation given in |Appendix A 



However, for (23) there is just a single turning point and now also the subdominant 



potential Qi{x) given by (25) not discussed in [TT]. The way subdominant terms such 
as Q2 enter is via the hierarchy of transport equations that arises in the JWKB method. 
To get the first correction the second or higher transport equations must be considered. 
In our case, Q 2 is present at the second transport equation. 

In the typical case without boundary several results are known for subdominant 
corrections to the leading Hamiltonian. For instance [23] gives a systematic discussion 
of ^-corrections considering a general series expansion, the Weyl symbol of the 
corresponding Hamilton operator allowing for momentum terms as well. For a general 
reference on JWKB-corrections see [22] ■ In our case, however, we must also consider 
the presence of a boundary xq. 

The result we find in Appendix A is a second order correction of the action 



I = I + e 2 I 2 , (55) 

where formally I2 is found by integrating I' 2 from the second transport equation. Just as 
the leading amplitude of the JWKB-function determined by the first transport equation 
diverges near the turning point, also the second transport equation leads to divergences. 
In our case the subdominant potential Q2 remains well behaved there. The singularities 
can be dealt with by for example transforming the problem to the Airy equation as in 



TT] and is given in detail in Appendix A 



To state the final result for the phase I2 we shall need the expansion of the potential 
Q in the coordinate x around the turning point a: 

Q(x) « a(x-a) + f3 (x - a) 2 . (56) 

Thus in the transition region x ~ a there is an approximate solution having the 
functional form of an Airy function which in the region x < a takes the form of a 
sine function. That region also possesses a standard oscillatory JWKB solution 

w osc (x) ~ Im [Cexp(il)] (57) 
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with the action / calculated up to a point before the turning point (a — fi with \i — > 0+) 



a— 11 



I(x) = / I'{t)dt. (5* 



We next match this oscillatory solution to the regular Airy-solution by adjusting the 
pre-factor C. This matching modifies the second order phase to 



lim ( 5 ^°) + 1 n Q " {t) ^ 

/3 a- 3 / 2 1 f a Q 2 (t) 



(59) 



12^ 2^ o (-g(t))va 

The first three terms are stated in [TTJ whereas the final term is the correction from the 



subdominant potential Q 2 . In (59) the first two terms and the last correspond to the 



result obtained by formally integrating the second transport equation for I 2 and doing a 



partial integration, see Appendix A . In this context, the second term diverges as Q(x) 



has a zero at the turning point. Given the expansion of Q(x) by (56), the third term 
coming from the matching exactly cancels this divergence and can be thought of as a 
regularization. 



Using computer algebra and further manipulations (59) reduces in the present case 

to 

/ 2 = Xo(3a4 + 2(x °- 3) : / 2 + a: ° ) + ^-latanf / X ° ) . (60) 
24a 2 (a 2 - x 2 f 2 16 8 \y/^~^J 1 > 

For an analytical check of the phase in the case of Xq = see |Appendix A 

As a numerical example, we continue that of the end of section [4] with 
I = 99.6428945787050, m = 70 and 8q = tt/3. Including the correction I 2 for the action 
leads to a scattering phase in units of 27r: Q m /(2ir) = 3.999993. Further numerical 
checks are given in Appendix A. 3. 3 



7.2. Effect of correction on the the trace formula 

The effect of the subdominant potential is to alter the saddle point position itlq to 

m* = mo + e 2 m 2 (61) 
to second order but this only affects the phases to fourth order: put 

S = 2n{I + e 2 I 2 ) = S + e 2 S 2 . (62) 
Then as the original phase function is stationary 

S* = S (m ) + (d m S (m ))e 2 m 2 + 0(e 4 ) = S* + 0(e 4 ) , (63) 
whereas the new part leads to 

S* 2 =e 2 S 2 (m ). (64) 

Finally, there is a contribution coming from a saddle point correction of the original 
integral, i.e. by expanding to fourth order in m around m , S'g 4 '*. There also other 
fourth order moments. However, these are of lower order in e and can be neglected. 
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At the point of stationary phase 



dm exp(— inO r 



2vre 

I Q"\ 
Po I 



dm exp 



S 



1/2 



exp I in 



7T 



7T 



exp 



So 
e 



eS 2 



exp 



S, 



(4) 



8S 



(2)2 



.(65) 



For a given orbit family, the factor 5*2 only alters the phase whereas the latter factor 
also alters both modulus and phase. 

The density of states is calculated using (33): essentially a derivative with respect 
to k = 1/e is performed on (65). Furthermore the derivatives of the action coefficients 
Si are pertaining to m and for these 



d 2 m d 

dm a da 

In particular Sq ~ e 2 . Thus apart from factors 

Sq 



Tr § ' 



~ e exp 



S 2 e + £e 



with 



and 



5, 



(4) 



a>i = «2 



v/2 8 ^ 2)2 ' 



Hence the perturbative correction to the density of states goes as 
~ <9 K (Tr §*") ~ i/^'S'oexp cos ^•S'o + 

where 

AS = S2 — ax 
Incorporating all pre-factors and constants yields 



1 

2^ 



p{k) 



00 00 



EE 

M=l n=l 



dA(p 



exp 



cos I nxK 



7T 7T AS\ 
U n + 7 + 

2 4 k y 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 



From (62) we notice the dependence on the number of repeats in the Si ~ n. Therefore 
AS = S2 + 0(n _1 ), so S2 = 2nl 2 with J 2 given by (60) in the case of the spherical cap 
remains the most important for high repeats. 



8. Summary and discussion 

We derived a trace formula for the spectral fluctuations of a spherical cap of general 
opening angle. The fluctuations in the spectral density were found to be governed 
by geodesic polygons. The method used is based on scattering referred to as inside- 
outside duality with a discussion of the explicit scattering states. Next the exact case 
of a hemisphere was introduced and there the spectral fluctuations were found to be 
governed by a single orbit of stronger weight than in the non-hemisphere case. The 
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leading result for general opening angles was compared to the hemispherical case at 
relatively high wave numbers kR « 45. Although non-uniform in nature, the general 
trace formula without JWKB-correction in the limit of opening angles close to 7r/2 also 
exhibits peaks at positions corresponding to the hemispherical cap's density of states. 
Finally, the leading JWKB-correction was discussed and incorporated. 

We used inside-outside duality in the derivation of the trace formula. The scattering 
states needed for this are perhaps not so familiar to physicists but these functions 
have already been introduced in seismology [9]. Furthermore inside-outside duality 
for billiards in curved spaces has been discussed recently in the thesis [12]. Despite 
that inside-outside duality is just one possible method for this problem of a spherical 
cap, it should also be mentioned that the method allows for generalizations and 
clear interpretations in more complicated cases such as systems of partial differential 
equations. 

Thus, although the leading result for the spectral density in principle can be 
obtained by a general theory of symmetry reduced trace formulae by [T5l [16] it is not 
always obvious how to generalize these results in applications for wave equations different 
from the Schrodinger equation: which reflection coefficients to use and what is the proper 
concept of the anholonomy entering for the overall amplitude of a family of orbits in 
the presence of for example ray-splitting. By construction, the scattering formalism 
automatically yields unitary reflection coefficients related to for example probability or 
energy flux. Furthermore, for the present work we have used the scattering method 
as a vehicle to go beyond the leading results and incorporate JWKB-corrections. In 
particular, such corrections to the spectral density cannot at present be obtained from 
the works of [TBI [TB] which pertain to the leading result only. We should also mention, 
that the scattering method can be extended to general shapes without symmetries by 
considering the scattering of suitable exterior states by for example attaching wave- 
guiding leads [3]. 

We expect the scattering method on spheres as discussed in this article and [T2] 
could be generalized to the case of a sphere with multiple circular holes. Here, addition 
formulae for the Legendre polynomials for non-integer angular momentum I [26] would 
lead to multiple-scattering expansions similar to those in flat space between discs or 
spheres [27]. 

As mentioned the work in this article has partly been motivated by Ellegaard's 
experiments on elastic shell caps. His group has studied plates and three-dimensional 
elastic resonators see [28, 29, 30J but is at the present time of writing investigating 
shells as well. In this context, the work in this article on the curved scalar Helmholtz 
equation and those in the flat case [UJ [311 E2] for two-dimensional elasticity show that 
derivations of trace formulae in more general settings are possible. Thus the combined 
case of curved elasticity [33], i.e. elastic shell caps, is indeed one generalization. 
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Appendix A. The JWKB solution 

Appendix A.l. Hierarchy of JWKB equations 
Using the oscillatory ansatz 

u> = Cexp {% (^- + I 1 + I 2 e) ) (A.l) 



e 

for the second order ordinary differential equation in the notation of [TT] 

e V = (Q + Q 2 e 2 )w (A.2) 

gives the eikonal equation 

(Q 2 + Q = (A.3) 

and the transport equation: 

- 2I[I' + < = (A.4) 
respective the second transport equation 

- 2I' 2 I' Q + ill ~ {1'if ~ Q2 = ■ (A.5) 
containing the correction from the subdominant potential. 

Appendix A.2. Leading order 

We fix the action Iq by integrating between x and the turning point a with the result 



(28). Likewise I x = zlog(-Q) 1/4 . From (A.l) and (22), a real JWKB solution in the 



oscillatory region away from the turning point is given up to two unknown constants A 
and 5 as 

u = A (1 - x 2 ) 1 / 2 (-Q)- 1/4 cos ^ + ^ = A (a 2 - x 2 )- 1 ^ cos + ^ (A.6) 

where A is an overall amplitude and 5 is a phase shift. 

One way of fixing A and 5 is by matching the wave function to an exact solution, 
here taken as 7rY; m (#, = 0) with 



me, * = o) = (-D™ ^a±lJ^±a , (a.7) 



at the point x = [H |26] using: 



2 m Y( l+m+1 ) / 77- 
^(0) = ^ f ^jcos(-(/ + m)) (Ai 
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and 



2 m+l jy H-m+2 -v /?r \ 

p r'(0) = F | ] _|_Z sin (_ (/ + m) ) . (A.9) 

To proceed with the match, the gamma functions in the normalization of the spherical 
harmonics are expressed via the duplication formula [BJ: 

T(2z) = -±= 2 2 - 1 / 2 r(z)T (z + ^ (A.10) 



and subsequently ratios of gamma functions are approximated with 

z a ~ b (A.11) 



T(z + a) 



T(z + b) 

valid for large z, see [BJ. As in the main text I and m are assumed asymptotically large 
with a fixed ratio. 

Next, the leading semiclassical phase at x = evaluates to 

/o = |(l-cosV) (A.12) 

which together with l/e~Z + l/2 enters for the shift 5. Finally cosip/e = m to all 
orders by definition (21). 

Thus, for A = 1 and 5 = — 7r/4 the wave function matches 

M^7rF im (M = 0). (A. 13) 

In particular, the condition for an eigenmode is 

= -+n7T (A.14) 

e 4 2 v ; 

for rz e Z. 



Appendix A. 3. The first correction 

We turn to the case where the sub-leading phase J 2 is included. We distinguish between 
the oscillatory region (osc), the transition region (trans) and the decaying region (exp). 
We shall not be concerned with the decaying region discussed in detail in [TJ] (chapter 
10.7). Instead, we focus on the oscillatory case from which it will be clear that the 
subdominant potential Q2 can be included in the discussion as well. 



Appendix A. 3.1. Phase We first remark, that we can always relate the solution at the 
boundary point x with that at x close to the turning point. The connection is by 
multiplication of the JWKB-factor 

uw(£o) = Im [exj)(iAI(x ,x)) w osc (x)] . (A. 15) 

with 

f x dl 

AI(x,x Q ) = ds — . (A.16) 
Jx ds 



Thus eventually all actions are continued further to the boundary as in (59). 
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Hence, the transport equations can be integrated formally, at least outside the 
transition region, whereas these diverge at the turning point and require discussion. In 
detail, to get the action to first order in e we need to integrate 

P _ __ Q2 Q" , 5 Q' 2 



2(-Q)V2 &(-Q)W 32 (-Q) 5 / 2 

* +^UAx*\ (A.17) 



2(-Q)V2 48(-Q) 3 / 2 48 dx\{-Q) 3 / 



as follows from (A. 5) with I' = +(—Q) 1 ^ 2 and I[ = z(log(— Q))' /A after some calculation. 
But, when this is integrated from x up to a — fi (with \x —>■ 0+) the last two terms diverge 
algebraically with respect to \i. 

We turn to the region near the turning point. The oscillatory solution is written in 
complex form as 

(x) = Cexp (i J " dsl'(s)j (A.18) 

We shall expand in the distance from the turning point 

Ax = x-a (A.19) 



with Aa; < corresponding to the oscillatory region. Thus from (56) 

-Q = \Ax\ all- -\Ax\) . (A.20) 



a 

The leading action 1$ is found from 



(-Q) 1 / 2 ~ \Ax\ 1/2 a^ 2 (l - |^|Ax|) (A.21) 

so 

The next action I\ is imaginary and contributes to the amplitude of the JWKB-function 
with the factor 



Q)- 1 / 4 ~ |Ax|" 1/4 a" 1 / 4 + ~|Ax|^ . (A.23) 



The second order action J 2 follows from integration of (A.17). First, 

Q" P IA _,_3/a 



Ax\~ 6/Z (A.24) 



48(-Q) 3 / 2 24a 3 / 2 
so 

Second, the exact term given by a derivative becomes 
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Finally the subdominant potential term with Q 2 = 7 + 0(Ax) (not having a zero at the 
turning point) becomes 

Q2 7 1 A _i-i/2 



Ax\~ 1/Z (A.27) 



2(-Q)V2 2aV2 
and leads to the regular 

"^(-^)^(i^ 1/2 -" I/2 '- < A - 28 > 

On the other hand, in the transition region x ~ a with 2 < a there is also a 
solution Wtrans in the form of the Airy function formula (10.7.7)). The Airy form 

is recognized in (A. 2) by the transformation 

x _ a = e 2/3 a -l/3 t _ (A>29) 

Then a perturbed Airy equation arises 

= (t + e^a-^^)w trans . (A.30) 

The leading perturbation due to Q is included whereas that from Q 2 does not play a 
role at this order. Consequently 

K'-SM^hS)- (A - 31) 

We now expand this Airy function in the region x < a to second order in an asymptotic 
series of oscillatory form: 

Ai(-z) ~ sin (C + I) - cos (C + I) - sin ( C + \ - ^ /2 ) (A.32) 

for large 2 with 

= n4/V 2 / 3 IA.d f 1 - ^ 

5a 



: = a^e-^IAarl ( 1 - ^ ) (A.33) 



and 



C = ? ^3/2 „ 2 al/2e _ 1|Ax|3 /2 ( x _ 3£jA*[\ . (A.34) 
3 3 \ 10a / 



Calculations similar as those for w osc give 

Wtrans(x) ~ 



Im 



5/2' 



L>VV 1 + e^cT^IAzr^exp i -a^V 1 \Ax\ a/z - 



4a 11 \ \ 3 \ 1 1 10a 



" 1 a .1-3/2-1/2. ^ , 3/3|Aa;|^ 



+ --— \Axr ,z a- i/z e 1 + 



, , , . (A.35) 
4 48' 1 V 10a y' ' 

This result for w tra ns is matched with the previous oscillatory form w osc . First, we 

1 — 1/2 

note that provided we neglect terms of order e|Ax| there is a complete match in 
functional form with respect to \ Ax\ in the oscillatory region close to the turning point. 
In particular, Q 2 does not play a role in the matching at this order. Second, as w tra ns 
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is independent of /i all dependence on /i must be removed in w osc . This is done by 
adjusting the pre-factor C of the standard JWKB-solution w osc = Im(C exp(if)). In this 
way singular terms in /i are absorbed. Thus for example also the last term's singularity 
in (A. 17) proportional to fx^ 3 ^ 2 is removed. 

Still in our case, the contribution from the subdominant Q2 to the second correction 
to the phase is regular and eventually gives a correction to the result stated in [11] for 
the full phase (i.e. all the way to the boundary x ): 



Ah 



a Q2(s) 



ds 



Q 2 (s) 



-Q(s))V- 



-ds . 



(A.36) 



Likewise the integration of the remaining terms in (A. 17) are seen to correspond to 
those given in (59). 

By construction the phase is now correct to first order in e. This can be verified at 

(A.37) 



the point x = as follows: There the phase I2 given by (60) evaluates to 

h = ~. 
16 

However, using 



i- l --r 2 +o{r 3 ) 



and 1/e = I + \ - \ T 1 + ±- V 2 + 0(r 3 ) (A.38) 
2 o 16 



and (|A.12|) as well as (|2l|) on the oscillatory part 

'7T 



if} oc COS 



cos 



(/-m) + 0(r 3 )) 



(A.39) 

<xY lm (0 = n/2,0) + O(r 3 ) 
shows the agreement in the phase even to second order in For further numerical 



checks see |Appendix A. 3. 3 



Appendix A. 3. 2. Amplitude Even though our main interest is a JWKB-resonance 
condition which only depends on the phase we note that similar calculations are possible 
for the amplitude function. These show agreement between the JWKB result and the 
exact to first order 



Appendix A. 3. 3. JWKB-phase at eigenmodes The table Al (opening angle 8q = 7r/3) 
shows the deviation of the JWKB-scattering phase measured in units of 2tc from the 
nearest integer (denoted [0/27r]) normalized with that integer, i.e. 

0/2tt - [0/2tt] 



[0/27T] ■ 

Apart from the value at (m, I) 



(A.40) 



^13,31.74) an overall improvement is observed. 
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Table Al. Scattering phase in units of 2n: the relative deviation from integer (A. 40 1 
at 8q — 7r/3. 



m 



I 



Leading Corrected 



m I 



Leading Corrected 







1 
1 
1 

2 
2 
2 
3 
3 
3 
4 
4 
5 
5 
6 
6 
7 
7 



1.777 
19.75 
37.75 
3.196 
21.24 
39.24 
7.622 
25.71 
43.73 
8.991 
27.16 
45.20 
19.53 
37.64 
17.87 
36.06 
16.14 
34.47 
20.59 
38.91 



-9.43 x lO" 3 
-1.32 x 10~ 4 
-3.76 x 10~ 5 
-3.79 x 10~ 3 
-1.22 x 10~ 4 
-3.61 x 10~ 5 
-8.31 x 10~ 4 
-8.69 x 10~ 5 
-2.99 x 10~ 5 
-4.80 x 10~ 4 
-8.00 x lO" 5 
-2.87 x 10~ 5 
-1.51 x 10~ 4 
-4.23 x 10~ 5 
-1.58 x 10~ 4 
-4.67 x lO" 5 
-9.84 x 10~ 5 
-5.11 x 10~ 5 
-8.76 x 10~ 5 
-4.03 x 10" 5 



-1.01 x 10 
-2.93 x 10 
-2.36 x 10 
-6.68 x 10 
-2.45 x 10 
-2.12 x 10 
-3.39 x 10 
-1.30 x 10 
-1.46 x 10 
-3.35 x 10 
-1.21 x 10 
-1.37 x 10 
-7.47 x 10 
-3.33 x 10 
-1.83 x 10 
-4.65 x 10 
-6.02 x 10 
-6.86 x 10 
-1.81 x 10 
-4.44 x 10 



-3 
-7 

-8 
-4 
-7 
-8 
-5 
-7 
-8 
-5 
-7 
-8 
-7 
-8 
-6 
-8 
-6 
-8 
-6 
-8 



9 

10 
10 
11 

12 
12 
13 
13 
15 
16 
17 
18 
19 
20 
21 
23 
25 
28 



21.94 
40.31 
26.40 
18.14 
37.00 
29.08 
20.68 
39.74 
31.74 
34.88 
24.46 
21.99 
23.21 
24.42 
29.45 
34.22 
38.87 
34.40 
36.86 
36.46 



-1 



-5.77 x 10 
-3.73 x 10 
-5.28 x 10 
6.75 x 10 
-3.95 x 10~ 5 
-2.42 x lO" 5 
8.30 x 10~ 4 
-3.06 x 10~ 5 
4.69 x 10~ 8 
-1.89 x 10 
1.01 x 10 
6.61 x 10 
6.72 x 10 
6.82 x 10 
1.17 x 10 
3.81 x 10 
1.59 x 10 
1.29 x 10 
1.34 x 10 
7.48 x 10 



-5 
-3 
-3 
-3 
-3 
3 
-1 
-1 
3 
-3 
3 



-1.81 x 10" 6 
-4.37 x lO" 8 
-7.20 x 10~ 7 
-3.30 x 10~ 5 
-1.05 x 10~ 7 
-7.19 x 10~ 7 
-3.29 x 10~ 5 
-1.05 x 10~ 7 
-7.18 x 10~ 7 
-3.39 x 10~ 7 
-3.27 x lO" 5 
-6.23 x 10~ 4 
-6.22 x 10~ 4 
-6.21 x 10~ 4 
-3.25 x 10~ 5 
-5.93 x 10- 6 
-1.80 x 10~ 6 
-3.23 x 10~ 5 
-3.22 x 10~ 5 
-6.14 x 10~ 4 



Appendix B. Scattering states 

Appendix B.l. Exact scattering states 

One possible set of exact scattering states for the sphere S 2 are those used extensively 
in geophysics under the name traveling waves in the discussion of the Green's function 
and associated surface waves [9]. We define them normalized: 



# )m = J r(l + m + l) [ Pr{x) ^ffii*) ) exp(im0) . (B.l] 



The opposite sign in (B.l) will be explained in the following. 



With this choice each state has constant total flux F through circles of constant 6: 
J • ds = I J § ds (B.2) 



e=cst 



Families of spherical caps: spectra and ray limit 



27 



with als = R sin 9 dtp. The corresponding probability current is proportional to (ignoring 
factors of H and mass): 
1 



Ji 



2i 



(+)m 



C.C. 



The directional derivative is rewritten as 
Id sin 9 d 

On 



Rd9 R dx ' 
where x = cos 9. The current in terms of the Wronskian becomes 
2 r(/-m + l) 



Js 



nR r(l + m+l] 



Wr(Pr,QT) sin9 



and hence the flux 



F = 4 sin 2 9 Wr(P™, Q" 



T(l-m + 1) 



(B.3) 
(B.4) 
(B.5) 

(B.6) 



T(/ + m + 1) ' 

Using the Wronskian [26] (p. 145: formula (25) in chapter 3.4 for Legendre functions on 
the cut: — 1 < x < 1 ): 



Wr(P-(x),QrW) 



(i - x 2 )r 



' Z— m+2 \-p/ ?— m+l ' 



(B.7) 



with the duplication formula for the Gamma function (A. 10) for the numerator and 



denominator in (B.7) then gives 
F = 4 > 



(B.8) 



corresponding to an outgoing state. The constancy of the flux reflects probability 
conservation. 



Appendix B.2. Semiclassical scattering states 



As in Appendix A we match the JWKB traveling states to the exact scattering states 
with the result 

7 tT 



{a 2_ x 2 r i/A 



exp ±z 



4 



using (A.8) and 



l+m+l- 



Qr(o) = -2™- i ^ 2 A_y sin (^( Z+m 



(B.9) 



(B.10) 



Hence, the scattering states lead to Dirichlet scattering phases given asymptotically by 



(-)m 



in 



(+)m 



CXp \ - I I ~~ + f ) ) = " '^Pf-'^m) ^ ' S '» 



(B.ll) 



with O m given by (30) to leading order 
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